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The production of rf mesons in photon- and hadron-induced reactions has been revisited in 
view of the recent additions of high-precision data to the world data base. Based on an effective 
Lagrangian approach, we have performed a combined analysis of the free and quasi-free jN ~ > rj'N, 
NN — > NNrj', and ttN — > rj'N reactions. Considering spin- 1/2 and -3/2 resonances, we found that 
a set of above-threshold resonances {Su, Pn, P13}, with fitted mass values of about Mr = 1925, 
2130, and 2050 MeV, respectively, and the four-star sub-threshold Pi3(1720) resonance reproduce 
best all existing data for the rf production processes in the resonance-energy region considered in 
this work. All three above-threshold resonances found in the present analysis are essential and 
indispensable for the good quality of the present fits. 

PACS numbers: 25.20. Lj, 13.60.Le, 13.75.-n, 14.20.Gk 
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I. INTRODUCTION 

A wealth of interesting physics can be obtained from 
studying the production processes involving the rf me- 
son one of the primary interests being here that 
such production processes may help one extract infor- 
mation on nucleon resonances that cannot be obtained 
from pion reactions. In fact, current knowledge of most 
of the nucleon resonances is mainly due to the study of 
irN scattering and/or pion photoproduction off the nu- 
cleon. Since the rf meson is much heavier than the pion, 
rf meson-production processes near threshold, therefore, 
are well suited for investigating high-mass resonances 
in low partial-wave states. Furthermore, reaction pro- 
cesses such as rj' photoproduction provide opportunities 
to study, in particular, those resonances that couple only 
weakly to pions. This may help in providing a better 
understanding of the so-called "missing resonances" pre- 
dicted by quark models, but not found in more traditional 
pion-production reactions Q. 

In view of the relatively low production rate, until re- 
cently there existed only a very limited number of ex- 
perimental studies of rf production reactions. This lim- 
ited experimental information was reflected in the rel- 
atively low number of theoretical investigations of such 
reactions. For an account of the pre-2004 situation, see 
Ref. pj. However, the situation has changed in the past 
few years, especially in rf photoproduction, where high- 
precision data for both nucleon and deuteron targets have 
become available Also, the pp and prf invariant- 

mass distribution data in the pp — > pprf reaction are now 
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available Ml , in addition to the cross-section data 
Upper limits for the total cross sections in pn — > pnrf 
have also been reported [ill ] . Given the present situation, 
with much higher-quality data than were available in the 
past, we revisit here the production of rf and we perform 
a combined analysis of the reaction channels 7 A?" — > rj'N, 
NN -> NNrf, and irN -> rj'N. 

To analyze the new photoproduction data [4}46|, we 
use here the same approach as employed in Ref. [12| for 
the analysis of the earlier CLAS data Q. In particular, 
we include spin-1/2 and -3/2 resonances with parameters 
determined from best fits to the data. Since the present 
approach is phcnomcnological, our strategy is to consider 
the minimum number of resonances that allows us to re- 
produce the available data with acceptable fit accuracy. 
This is quite different in spirit from a recent analysis of 
the new data [4-6] by Zhong and Zhao [l3| who work 
within a quark-model approach and consider all possible 
nucleon resonances up to the n = 3 harmonic-oscillator 
shell. For the quasi- free photoproduction processes [6j, 
we account for the Fermi motion of the nucleon by fold- 
ing the cross section of the free process with the momen- 
tum distribution of the nucleon inside the deuteron. The 
analysis of the reaction NN — > NNrf is done follow- 
ing Ref. [l| within a distorted-wave Born approximation 
(DWBA) in which both the initial- and final-state NN 
interaction is taken into account explicitly. 

This paper is organized as follows: In Sec. UTJ we de- 
scribe our overall strategy for performing our combined 
analysis of the photon- and hadron-induced reactions and 
we also provide some general remarks concerning the res- 
onances required to reproduce the available data. In 
Sec. IIII1 the present results for the free ^p — > rfp re- 
action are discussed in conjunction with the most recent 
high-precision CLAS and CBELSA/TAPS [| data. 
In Sec. IIV1 the analysis of the quasi-free 7^ — > rfp and 



2 



-fn —> rjn reactions is presented. Section [V] contains 
the present results for the NN — > NNrj' reaction and 
in Sec. IVII we present the results for the pion-induced 
i]' production TtN — > rj'N. Our summarizing assess- 
ment is given in Sec. IVIII In Appendix [XJ our model 
for photoproduction as well as for the hadron-induced 
reactions NN —> NNn' and irN —> n'N are described 
briefly for completeness. Some details of how the Fermi 
motion of the nucleon inside the deuteron is taken into 
account in the present work for describing the quasi-free 
photoproduction processes are also given in Appendix [A"l 
Appendix [Bl contains the Lagrangians, form factors, and 
propagators that define the individual amplitudes in our 
model. 



II. GENERAL PROCEDURE AND FINDINGS 

In the following sections, we present separate discus- 
sions of our results for photon- and for hadron-induced 
reactions because this allows us to keep the discussions 
focused on the reactions at hand. We emphasize, how- 
ever, that the results are based on simultaneous fits to 
the available data for the reaction processes considered. 

The model assumptions used in describing these re- 
actions consistently with each other are given in Ap- 
pendix^ In addition to resonance-current contributions, 
we consider the nucleonic and meson-exchange (p and 
lj) currents. For resonances, in particular, the strat- 
egy adopted here is to introduce as few resonances as 
possible, with parameters adjusted to simultaneously re- 
produce the data for the 7N rj'N, NN NNrf, 
and irN — > vj ' N reactions via a least-square minimization 
procedure. In the present work, we restrict ourselves to 
spin-1/2 and -3/2 resonances. 

The simultaneous treatment of the available reaction 
data provides insights that cannot be obtained by fit- 
ting each reaction and/or each data set individually. For 
example, considering only the free *yp — > rj'p photopro- 
duction data from the CBELS A/TAPS experiments 
we can obtain good fits (with /N ~ 1) with only 
one resonance. By contrast, a similarly good fit to the 
corresponding high-precision CLAS data [U requires at 
least three resonances. Considering two above-threshold 
resonances instead, one obtains x 2 /N = 1.7 at best. 
Moreover, we find that any one of the above-threshold 
three-resonance sets {Pn, -P13, -D13}, {Sii, Pn, P13}, and 
{S11, P13, -D13} can reproduce the data equally well. The 
remaining combination {5*11, Pn, ^13} is ruled out by the 
data according to our good-fit criterion thus indicating 
the need for the P13 resonance for reproducing the data. 
These findings indicate that in 7^ — > rfp the cross-section 
data alone cannot constrain the set of resonances, even 
with high-precision data. By contrast, the three accept- 
able sets of resonances given above yield quite distinct re- 
sults for the spin observables, in particular, for the beam 
asymmetry. These findings corroborate the conclusions 
of our earlier work [121. 



The ambiguity just discussed with respect to the three 
resonance sets found acceptable in the free photopro- 
duction process "fp — » rj'p is completely removed once 
we include other reactions in our analysis. The set 
{Pn, P13, P13} is ruled out by the quasi-free — > r/'n 
reaction and the set {S11, P13, P13} by the JVJV ->■ NN-q' 
reaction. The only remaining acceptable set, therefore, is 
{S11, Pn, P13}. In addition, we find that to obtain a good 
description of NN — > NNn', in particular, we need an 
additional &ek>w-threshold resonance whose inclusion has 
no bearing on the quality of any of the results for other 
reactions considered in this work. In summary, there- 
fore, a good overall description of all reaction processes 
considered in this work is obtained with a (minimum) 
set of three above-threshold resonances, {Sn, Pn, P13}, 
with the masses of about Mr = 1925, 2130, and 2050 
MeV, respectively, and one below-threshold Pi 3 reso- 
nance. The latter is the four-star Pi3(1720) quoted in 
PDG 14]. The masses of these above-threshold reso- 
nances are very well constrained by the existing data, in 
particular, by the CLAS photoproduction data [||. The 
CLAS data also constrain reasonably the total widths of 
these resonances, although to a lesser extent. Hereafter, 
we shall use the notation 5 U (1925), Pu(2130), Pi 3 (2050) 
to identify these above-threshold resonances, even though 
the fitted mass values for the various scenarios discussed 
below (see Table HJ may differ slightly from these val- 
ues. The resonances found in our analysis may be ten- 
tatively identified with the corresponding 511(1895)**, 
Pn(2100)*, and Pi 3 (2040)* resonances listed in PDG 
[l<lj . We will discuss the constraints imposed by the in- 
dividual reactions in more detail in subsequent sections. 

III. FREE 7p rj'p 

In this section, we present our results for the free 
7P — > r/'p reaction. First, we address the issue of possible 
discrepancies in the data sets from the CLAS [3] and the 
CBELSA/TAPS @ collaborations. Then, we discuss the 
present analysis of these data. 



A. Comparison of the CLAS and CBELSA/TAPS 

data 

The most recent CLAS data [i[ and the 
CBELSA/TAPS data @ are compared in Fig. [Q 
As one can see, in general the new CLAS and the 
CBELSA/TAPS data are consistent for invariant en- 
ergies below ^fs ~ 2.0 GeV within their uncertainties. 
For higher energies, however, one sees considerable 
discrepancies, by factors as large as 3, between these two 
sets of data for 7/ emission angles away from the forward 
angles. We note that in addition to the statistical errors, 
the plots of Fig. [1] include the estimated systematic 
errors as quoted in Refs. 0, [f|. These systematic errors 
were not included in the data plots given in these 
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FIG. 1. (Color online) Comparison of the CLAS @] (blue solid circle) and CBELSA/TAPS Q (red open square) differential 
cross-section data for the free 7p — > r/'p reaction as a function of cos(#) (where 8 is the rf emission angle in the center-of- 
momentum frame) for invariant energies up to \fs = 2.35 GeV. The earlier CLAS data || (black solid square) are also shown. 
The numbers in parentheses denote the photon laboratory incident energy (left number) and the total energy (right number) of 
the system, with the upper number pair pertaining to the newer CLAS data, the second one to the CBELSA/TAPS experiment, 
and (where present) the lower one to the older CLAS data. Note that the energies for different experiments shown in the same 
panel are within ±10 MeV. 



references. The error bars in Fig. Q] were obtained 
by adding the systematic errors to the corresponding 
statistical errors in quadrature. For the present purpose 
of revealing the inconsistencies between these two data 
sets, this procedure may be justified, even though one 
might need a more thorough error analysis for a fully 
quantitative estimate of the total uncertainty in the 
data. We mention that for 77 photoproduction, an even 
more pronounced discrepancy between the CLAS and 



CBELSA/TAPS data was pointed out in Refs. 0, 
however, with no clear identification of the source of 
the discrepancy (in this respect, see also Ref. (l5j|). For 
comparison, we show in Fig. Q] also the earlier CLAS 
data Q which are seen to be consistent with the newer 
CLAS data. 

Leaving out the older CLAS data [3j , we show in Fig. [5] 
the same data as in Fig. [I] however, now with curves that 
result from our fit procedure to the new CLAS Q and 
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FIG. 2. (Color online) Same differential cross-section data as Fig. [T] (but without older CLAS data fjj), now with curves 
resulting from our fits discussed in Sec. IIII Bl The CLAS data and fit curves are shown as blue solid dots and blue curves, 
respectively, and the CBELS A/TAPS data and fit curves are depicted as red open squares and red curves, respectively. The 
curves are the results of a combined fit of the photon- and hadron-induced reactions data and both include the set of resonances 
{Pi3(1720), 5n(1925),Pn(2130), Pi3(2050)} with parameters given in Table [I] All other parameter values that influence the 
photoproduction reaction directly are also given in Table U 



the CBELSA/TAPS ® data. We wil l discu ss these fit re- 
sults in detail in the subsequent Sec. IIII Bl Here, we note 
that if we integrate our fit results for the corresponding 
angular distributions in Fig. [21 we obtain the total cross 
sections shown in Fig. [3] which exhibit markedly different 
behavior for the two data sets. First, these results clearly 
reveal an energy-dependent relative normalization prob- 
lem which increases with increasing energy. We checked, 
of course, how much this finding is influenced by the un- 
certainties of the fit procedure and we found that other 



fits to the data of similar quality, even polynomial fits, 
have no effect on our conclusion regarding the normal- 
ization problem. Again, a similar finding was observed 
for 77 photoproduction data [3, [j| , as was pointed out in 
Ref. Second, the CBELSA/TAPS data exhibit a 

pronounced peak structure around 2.05 GeV, while the 
total cross-section resulting from the fit to the CLAS data 
shows only a relatively flat bump at this energy. The ori- 
gin of this pronounced peak structure can be traced back 
to the measured cross sections in the CBELSA/TAPS 
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FIG. 3. (Color online) Predictions for the total cross sections 
in (free) 7p — >• rj'p obtained by integrating the correspond- 
ing fit results of the CLAS [J] (dashed blue line) and the 
CBELS A/TAPS |[ (solid red line) angular distribution data 
as shown in Fig. fusing the parameters of Table [I] The data 
are from CBELS A/TAPS 5] obtained by simply integrating 
the corresponding differential cross section. They are not in- 
cluded in the fit. 



data bin around 2.052 GeV. We will discuss this struc- 
ture in more detail in the following subsection. 



B. Analysis of the CLAS and CBELS A /TAPS data 

The discrepancy between the CLAS 4] and 
CBELSA/TAPS 3 data discussed in the previous 
subsection makes it difficult to consider both data sets 
as a single set of data for a combined numerical analysis. 
Since, at present, we have no clear reason to discard one 
data set in favor of the other, we are forced, therefore, 
to consider them separately. In view of their differences, 
we can expect that the resonance parameters extracted 
from these data sets may be quite different from each 
other. One of the purposes of this subsection is to see 
how different they are. 

As mentioned, Fig. [5] shows the independent resulting 
fit curves for the CLAS [4] and for the CBELSA/TAPS 
5] data. Both data sets are reproduced with very good 
fit quality of \ 2 /N = 0.62 and 0.56, respectively. The 
corresponding model parameter values are displayed in 
the two columns subsumed under "free p" in Table fl] 
The uncertainties in the resulting parameters are esti- 
mates arising from the uncertainties (error bars) associ- 
ated with the fitted experimental data points. In addi- 
tion to the resonance-mass and total-width values, the 
table shows the corresponding reduced helicity ampli- 



tudes J PniY Aj , where (5n v ' denotes the branching ratio 
to the decay channel Nrf and Aj stands for the helic- 
ity amplitude with spin j. The mass and total width of 
the (four-star) Pi3(1720) resonance have been fixed at 
the respective centroid values quoted in PDG [14|. Also, 
the fixed radiative decay branching ratio of (3 pl = 0.10% 
is well within the range of [0.05-0.25]% quoted in PDG 
for this resonance. The analysis of r\ production pro- 
cesses in Ref. [lj| yielded a value of /3 P7 ~ 0.12%. We 
recall that this sub-threshold resonance is required in de- 
scribing the NN — > NNrj' reaction but not in photo- 
production. As such, its contribution is negligible here 
(cf. Fig. HI and, consequently, its parameters are not 
well constrained by the present photoproduction data. In 
particular, the corresponding reduced helicity amplitudes 
quoted in Table U are subject to much larger uncertain- 
ties than those indicated there, since the extracted very 
small branching ratio /3jvi/ for this resonance (cf. Ta- 
ble IIII[) depends more sensitively on the assumed mass 
distribution (see discussion in Sec. [Vj) than the corre- 
sponding values for the above-threshold resonances. The 
systematic uncertainties arising from such an assumption 
are not taken into account in the error numbers given in 
Table HI Following Refs. 0, [Hj], we also set the radia- 
tive decay branching ratios for other resonances to be 
j3 pl = 0.2%. Note that in the present tree-level calcu- 
lation, the results are rather insensitive to the branch- 
ing ratios, since they enter the model only through the 
total decay widths in the resonance propagator (see Ap- 
pendix IB 3|) , in addition to the fact that the results are 
sensitive only to the product of the coupling constants 
gRNr/'gRN-y- In principle, a simultaneous consideration 
of photon- and relevant hadron-induced reactions would 
enable us to disentangle the individual factors contribut- 
ing to this product of coupling constants. However, the 
relatively poor quality of the currently existing data for 
one of the relevant reactions, ttN —> rj'N (cf. Fig. [T51dis- 
cussed in Sec. IVII below) . does not allow this. We also 
note that the pseudoscalar-pseudovector mixing parame- 
ters A (see Appendix lB 2\ are not well constrained by the 
data, as the corresponding uncertainties indicate. The 
fitted mass values (Mr) of the corresponding resonances 
found from the CLAS and from the CBELSA/TAPS data 
are very close to each other. The same is true for all but 
one of the total widths Tr of these resonances at the res- 
onance energies W = Mr. A marked exception is the 
Pi 3 resonance for which the CBELSA/TAPS data yield 
a much narrower width (cf. Table [J), however, with an 
associated uncertainty that is very large. We will dis- 
cuss this issue in more detail below in connection with 
the total cross-section results shown in Fig. |4j The pa- 
rameter values for the Sn resonance are practically the 
same for both data sets. We emphasize that the CLAS 
photoproduction data [i| constrain the masses of the 
above-threshold resonances very well. In contrast, the 
CBELSA/TAPS data @ by themselves cannot constrain 
these resonance masses because of the overdetermined fit 
(recall that CBELSA/TAPS data [5] can be fitted only 
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TABLE I. Model parameter values that directly affect the photo-production reaction obtained in a combined analysis of the 
photon- and hadron-induced reactions. The values in the columns labeled CLAS and CBELSA/TAPS subsumed under "free p" 
correspond to the fit results for the CLAS and CBELSA/TAPS [f| free proton data, respectively. The last row corresponds 
to the fit results for the CBELSA/TAPS quasi- free proton data [f| discussed in Sec. IIVI Values in boldface were kept fixed 
during the fitting procedure. (The fixed mass values for the quasi-free calculation are obtained as averages of the corresponding 
CLAS and CBELSA/TAPS values; see text for further explanation.) For the definition of the parameters, see Appendices |A1 
and|U The resonance mass and total width, Mr and Vr, are both in units of MeV, while the reduced helicity amplitude, 
yffi^Aj, is in units of 10" 3 GeV" 1/2 . 



free p quasi-free p 





CLAS 


CBELSA/TAPS 




X 2 /N 


0.65 


0.53 


0.77 


Q AT \Tn f 


1.00 ±0.06 


1.17 ± 0.31 


1.00 ±0.24 




0.53 ±0.06 


0.44 ± 0.22 


0.64 ± 0.24 


A v [MeV] 


1183 ±5 


1244 ± 35 


1221 ± 28 


h 


3.89 ± 0.18 


5.37 ± 1.57 


4.27 ±0.89 


Pis (1720) 








Mr 


1720 


1720 


1720 


Tr 


200 


200 


200 




0.09 ±0.03 


0.09 ±0.06 


0.06 ±0.11 


\l 3 JVn' A-Z /O 
V " JV '/ o/ z 


-0.16 ±0.05 


-0.13 ±0.09 


-0.03 ±0.06 


Pi 3 (2050) 








Mr 


2050 ± 4 


2045 ± 7 


2048 


Tr 


140 ± 10 


52 ± 184 


51 ± 241 




-5.71 ±0.17 


-2.02 ±0.26 


-3.14 ±0.43 




9.89 ± 0.30 


7.31 ± 0.93 


5.75 ± 0.79 


5n(1925) 








Mr 


1924 ± 4 


1926 ± 10 


1925 


Tr 


112 ±7 


99 ±23 


145 ± 45 


A 


i oo +0 -°° 

1,uu -0.06 


1 00+ 000 

±,uu -0.98 


100 +o.oo 

-0.95 


\J PNri' ^4-1/2 


-11.84 ±0.41 


-11.07 ± 1.43 


-19.93 ± 1.56 


Pn(2130) 








Mr 


2129 ± 5 


2123 ± 23 


2126 


Tr 


205 ± 12 


246 ± 54 


170 ± 178 


A 


1 00+ 00 

± -" KJ -0.04 


100+ o.oo 


± - uu -0.95 




-11.34 ±0.62 


-18.80 ±0.90 


-7.45 ±0.94 



with one above-threshold resonance). 

The differences in the corresponding individual cur- 
rent contributions are better seen in Fig. |4j Comparing 
the resonance values found in our fits with those quoted 
in PDG [3, our Pn(2130), Su(1925), and Pi 3 (2050) 
may perhaps be identified with the PDG resonances 
Pii(2100)*, 5x1(1895)**, and Pi 3 (2040)*, respectively. 
As mentioned in Sec. UH the cross section data alone in 
7P — > rfp cannot constrain the set of resonances uniquely. 
We emphasize, however, that acceptable fits of the cross- 
section data for this reaction cannot be obtained without 
the above-threshold Pi3(2050) resonance. 

Figure [4] displays the individual resonance contribu- 
tions to the predicted total cross sections obtained by 
integrating the differential cross section results shown in 
Fig. OH for the CLAS (left panel) and CBELSA/TAPS 



(right panel) data. Both the nucleonic and mesonic cur- 
rents yield similar contributions for both data sets, even 
though they are visibly larger for the CBELSA/TAPS 
data than for the CLAS data. For resonances, we 
see significant differences in their relative contributions. 
The Pn resonance contribution is much stronger for 
CBELSA/TAPS than for CLAS, which is responsible 
for making the full total cross section larger for ener- 
gies above W ~ 2.1 GeV. For both data sets the respec- 
tive S\\ resonance contributions alone are responsible for 
the sharp rise of the full total cross section (red curves) 
near threshold, and their corresponding parameter sets 
are practically identical. By contrast, the sub-threshold 
Pi3(1720) resonance contribution is negligible. 

The clear peak structure in the full total cross exhib- 
ited by the CBELSA/TAPS data is produced here by the 
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FIG. 4. (Color online) Total cross sections with individual (nucleonic, mesonic, and resonance) current contributions. The left 
panel results from the fit to the CLAS data while the right panel pertains to the CBELSA/TAPS data. The corresponding 
parameters are given in Tabled] The peak structure at 2.05 GeV for the latter data is solely due to the narrow Pi3(2050) 
resonance. 
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FIG. 5. (Color online) Effects from p and w exchanges. The red solid line is the constructive sum of the p (green dashed) and 
oj (black dash-dotted) contributions, while the magenta dotted line is the destructive sum of p and uj. 



above-threshold Pi 3 resonance which is much narrower 
for the CBELSA/TAPS data than for CLAS. However, 
the uncertainty associated with its width is very large, 
as shown in Table U and this width, therefore, is not 
well constrained. We recall here that the total cross- 
section data were not fitted, so the peak structure found 
in the calculated results is a consequence of the fit re- 
sults of the angular distribution data as shown in Fig. 



The origin of this large uncertainty in the width can be 
traced back to the CBELSA/TAPS measured differential 
cross-sections at one particular energy only, namely at 
W = 2.052 GeV, which coincides with the Pi 3 resonance 
whose position is fixed very well within our fits with very 
small errors for both data sets. In view of the large er- 
ror for the P13 width for the CBELSA/TAPS data, there 
is no physical significance in finding the peak structure 
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FIG. 6. (Color online) Comparison of the quasi-free 7p — > r/'p differential cross-section data from CBELSA/TAPS [6( (open 
green squares) with the corresponding free data from CBELSA/TAPS [f| (solid red squares) and from CLAS [|[ (solid blue 
circles). 



shown in Fig. U In fact, manually changing the width 
from 52 to 140 MeV, similar to the width found for the 
CLAS data, largely smoothes out the peak resulting in 
an overall y 2 /N only a few percent worse than what is 
reported in Table |TJ 

The main conclusion regarding the discrepancy be- 
tween the CLAS Q and CBELSA/TAPS @ data as ex- 
hibited in Fig. [T] is that the larger cross-section yield of 
the CBELSA/TAPS data at higher energies results in a 
larger P±i resonance contribution compared to the CLAS 
data (compare, in particular, the corresponding reduced 
helicity amplitudes in Table fl} . This alone largely leads 
to the enhancement of the CBELSA/TAPS total cross 
section over the CLAS results seen in Fig. [3] for energies 
above W ~ 2.1 GeV. For future analyses, it is impor- 
tant to resolve this discrepancy in the findings of the 
CLAS and CBELSA/TAPS experiments if one is to ob- 
tain more definitive answers about these above-threshold 
resonances. 

Figure [5] illustrates the contributions of the p- and 
w-meson exchanges to the total mesonic current. De- 
pending on their relative sign, one obtains constructive 
(red solid curves) or destructive (magenta dotted curves) 
interference. The motivation for showing this detail 
here is that the MAID group [l?} employ p and ui cou- 



pling constants in their Regge contribution correspond- 
ing to a destructive interference very close to the result 
shown here for a destructive interference. By contrast, in 
our present calculations we have a constructive interfer- 
ence between the p- and w-meson contributions following 
Refs. Obviously, one gets different resonance pa- 

rameters whether the interference for this mesonic back- 
ground is constructive or destructive. As can be seen 
from Fig. [SJ in the former case, the total mesonic cur- 
rent (red solid curves) is relatively large at forward an- 
gles, while in the latter case (magenta dotted curves), it 
is much smaller over the entire angular range. In our 
model, the p and w coupling constants (including the re- 
spective signs) are determined from the measured decay 
widths and the SU(3) symmetry considerations in con- 
junction with the OZI rule (cf. Refs. [H, H2]) which leads 
to a constructive interference. We were unable to iden- 
tify how the signs of the p- and w-meson contributions 
were determined in Rcf. 11711. 



IV. QUASI-FREE -yN -> rj'N 

In this section, the quasi-free photoproduction pro- 
cesses 7p — ¥ rj'p and jn — > rj'n are discussed. They are 



9 




-1.0 -0.5 0.0 0.5 1.0 -0.5 0.0 0.5 1.0 -0.5 0.0 0.5 1.0 -0.5 0.0 0.5 1.0 -0.5 0.0 0.5 1.0 -0.5 0.0 0.5 1.0 

cos(e) 



FIG. 7. (Color online) Comparison of the quasi-free jp — > rj'p differential cross-section CBELSA/TAPS data of Fig. [6j with 
theoretical results obtained via Eq. (1A8[) . The data [6| are represented here by solid (blue) squares. The (cyan) dashed and the 
(olive) short-dashed curves are obtained by respectively folding the CBELSA/TAPS and the CLAS fits of Fig. [2] The latter 
two folded fits have no extra parameters. The (red) solid curves provide the fit results of the quasi-free data, also obtained by 
using Eq. (|A8[) . The parameters of this fit are given in the right-most column of Table [I] 



evaluated by folding the corresponding cross sections for 
the free processes with the momentum distribution of the 
target nucleon in the deuteron as described in Appendix 
S[cf. Eq. (|A8]) ]. 

We first compare in Fig. [5] the quasi-free 7p — > rj'p 
cross-section data obtained by the CBELSA/TAPS Col- 
laboration [6J with the corresponding free jp rj'p data 
from Refs. [J, H|. We note that, as for the free jp —> r)'p 
data, the uncertainties in the quasi-free data shown here 
include both the statistical and systematic errors added 
in quadrature. The latter errors were quoted in Ref. @ 
but were not included in the data shown in that refer- 
ence. Here, one sees that the quasi-free data practically 
coincide with the free data within their uncertainties for 
most of the angles and energies, although one sees some 
tendency for stronger angular dependence at a few ener- 
gies. However, this does not necessarily mean that the 
free and quasi-free data are compatible with each other 
as discussed below. It is also interesting to note that, 
overall, the quasi- free CBELSA/TAPS data are more 
in line with the free CLAS data Q| than with the free 
CBELSA/TAPS data 0], a feature that has been also 
pointed out in Ref. @. 

If the quasi-free hypothesis holds true, one would ex- 



pect that the Fermi motion of the target nucleon in- 
side the deuteron smears out the energy dependence of 
the free cross section. This then should affect more the 
cross section at low energies where one usually observes a 
strong energy dependence in the corresponding free cross 
section (cf. Fig. [3]). At higher energies, where the en- 
ergy dependence of the corresponding free cross section 
becomes weaker, the Fermi-motion smearing should have 
very little effect on the cross section. Therefore, based 
on the free total cross section results shown in Fig. [HJ it 
is conceivable to expect the Fermi motion to affect the 
cross section up to about W ~ 2 GeV. However, this is 
not what we observe in Fig. [51 where the quasi-free and 
free data still coincide with each other down to the low- 
est energy (W = 1.935 GeV) for which the free data are 
available. It would be interesting to also have data for 
the free process for lower energies closer to the threshold 
energy (at 1.896 GeV) to see the energy region where 
Fermi motion is relevant in this reaction. 

Figure [7] shows the same quasi-free jp — > rj'p differ- 
ential cross-section data as in Fig. [6j but now with our 
results. The cyan dashed curves are obtained by simply 
folding the fit results of the free CBELSA/TAPS data [| 
shown in Fig. [5] with the momentum distribution of the 
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FIG. 8. (Color online) Fit results for the differential cross section of the quasi-free — > rj'n reaction as a function of 003(6*) 
(where 9 is the r/' emission angle in the center-of-momentum frame) for energies up to y/s = 2.35 GeV. The data are from 
CBELS A/TAPS [|. 



proton inside the deuteron according to Eq. (|A8I) . and 
the olive short-dashed curves correspond to the analo- 
gous results obtained for the free CLAS data There 
are no extra fit parameters here; the parameters are all 
pre-determined by our fits of the free — > rj'p reaction. 
First, we see that the differences between the two sets of 
results obtained by folding are much smaller than what 
we observed in Fig. [5] for the free process at low ener- 
gies. Here, the Fermi motion may be smearing out the 
differences observed there. Comparing with the data, 
we see that the present model predictions are overall 
fairly reasonable considering the fact that there are no 
free parameters to fit. We reiterate here, however, that 
the quasi-free data exhibit a stronger angular dependence 
for some energies, a finding already pointed out in con- 
nection with the discussion of Fig. [BJ Also, the present 
predictions seem to exhibit a slight tendency to under- 
estimate the lower-energy quasi-free data. This feature 
will be clearer as discussed below. The prediction corre- 
sponding to the fit result of the free CBELSA/TAPS data 
@ (cyan dashed curves) tends to overpredict the data at 
some higher energies. Overall, the fits of the folded free 
CLAS [| and CBELSA/TAPS @ data have an increased 
X 1 /N of 1.4 and 2.5, respectively. 

The red solid curves in Fig. [7] correspond to the fit re- 



sults of the quasi- free data also obtained via Eq. (|A8I) . 
The resulting parameters for this fit are shown in the 
right-most column of Table U for a direct comparison 
with those resulting from the free data fit. Here, the 
resonance mass values were fixed to be the average of the 
fit results of the free CLAS and CBELSA/TAPS data, 
since they are well determined by these free data. As 
can be seen, the fitted values of the other parameters are 
close to the corresponding values obtained from the fit 
of the free data as expected for the quasi-free process. 
The only notable difference is in the reduced helicity am- 
plitudes, yj fiNri'Aj, for the spin- 1/2 resonances. For the 
Pn(2130) resonance, the value is much closer to the free 
CLAS result, which is understandable because at higher 
energies, the quasi-free data are much closer to the free 
CLAS data than to the free CBELSA/TAPS data. This 
might indicate a possible normalization problem in the 
free CBELSA/TAPS data at higher energies. For the 
Sii(1925) resonance, the reduced helicity amplitude is 
almost a factor of 2 larger than the corresponding values 
extracted from the free data. As can be seen in Fig. [3 
this is also easy to understand; the fit simply tries to 
enhance the cross section at lower energies, where the re- 
sults obtained by folding the free cross sections without 
any fit parameters tend to underestimate the data. 
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Summarizing the discrepancies between quasi-free data 
and the Fermi-folded free cross sections, we find that the 
folded free results tend to underestimate the quasi-free 
data at lower energies and for higher energies, we find 
that, if anything, the folded free results are above the 
quasi-free results for some energies, in particular for the 
CBELSA/TAPS results. There might be two possible 
causes for this energy-dependent difference between the 
theoretical folding procedure and the corresponding ex- 
perimental analysis: (i) Our prescription for accounting 
for the Fermi motion is not quite adequate [sec Appendix 
IA 21 in particular, Eq. (|A8I) ]. This prescription, however, 
works quite well for r\ photoproduction [l8|, [l9| where 
there is a much stronger energy dependence in the cross 
section close to threshold. To test how sensitive the re- 
sults are on the details of the folding procedure, we have 
also employed the alternative prescription of Eq. (|A7[) 
which treats the total energy available to the 7p — > rj'p 
subsystem differently from Eq. (|A8j) (for details, see Ap- 
pendix |A~21 . and we found no appreciable differences, (ii) 
The quasi-free data contain additional nuclear effects at 
lower energies which cannot be adequately described by 
the simple folding procedure. In fact, our calculation 
shows in Fig. [7] that the effect of Fermi folding can be 
seen for energies up to about 2 GeV, while the direct 
comparison of quasi-free and free data in Fig. [6] do not 
show this effect for energies down to the lowest energy 
of 1.934 GeV available for the free data. We note in 
this context that to force the theoretical folding results 
to agree with the quasi-free data at lower energies, we 
would need to drastically cut the higher-momentum part 
of the deuteron wave function. Further investigation are 
necessary to find out which of these two possible causes 
applies here — or perhaps even a combination of both. 

One purpose of fitting the quasi-free data here is to 
have a quasi-free proton result on the same footing as 
the quasi-free neutron result, since the latter has to be 
fitted to the corresponding quasi-free data to fix the reso- 
nance transition electromagnetic couplings. In this way, 
quantities such as the ratio of the neutron and proton 
branching ratios will be free of possible unwanted effects 
which may distort the results otherwise. 

In Fig. [8] we show the results for the quasi-free —> 
rj'n reaction. As mentioned in Sec. HH incorporating this 
reaction in our combined analysis rules out the set of 
above-threshold resonances {Pn, P13, -D13} which other- 
wise for the free jp — > rj'p reaction reproduces the data 
just as well as the set {Su, P%i, -P13}. Overall, we re- 
produce the quasi-free "/n — > rj'n data reasonably well, 
with x 2 /N = 0.82. The calculation requires adjusting 
free resonance-neutron-7 (Rwy) coupling constants that 
need to be determined through the combined fit of the 
photon- and hadron- induced reaction data. Table |TT] dis- 
plays the resulting values. For the Pi3(1720) resonance, 
we obtain /3„ 7 = 0.016% which is at the upper limit of 
the range of [0.0 - 0.016]% quoted in PDG [lj|. We note 
that, as in the free jp — ¥ rj'p case, here the parame- 
ters associated with this sub-threshold resonance are not 



TABLE II. Electromagnetic couplings extracted from the 
CBELSA/TAPS quasi-free neutron data fl in a global fit 
of the photon- and hadron-induced reactions data. The cor- 
responding branching ratios of resonances decaying into nj 
and P7, Pn-y/Pp-,, are listed in the last row. The sign of this 
ratio reflects the relative sign in the corresponding neutron 
and proton electromagnetic coupling constants. 





Pi 3 (1720) 


Pis (2050) 


5u(1925) 


Pn(2130) 


\J Pnti'A-l/2 


0.04 


0.94 


15.54 


7.60 




-0.00 


-1.64 








0.32 


-0.09 


-0.61 


-3.06 



well constrained by the data. In Table [TT1 we also give 
the ratio of the neutron-to-proton branching ratios. The 
sign of this ratio reflects the relative sign in the corre- 
sponding neutron and proton electromagnetic coupling 
constants. Although the branching ratios may be subject 
to a considerable ambiguity since, as discussed in the pre- 
vious section, only the product of the coupling constants 
9RNri' 9RN-y is well determined in the present calculation, 
the ratio /3„ 7 //3 P7 is free of such an ambiguity. 



V. NN NNrj' 

In this section, we discuss our results for the nucleon- 
induced reaction JVJV — > NNrj' obtained from the com- 
bined analysis of this reaction together with the irN — > 
rl N and 7^ — > rj' N reactions. As mentioned in Sec. HJ 
we find that among the spin-1/2 and -3/2 resonances 
considered in this work, the set of above-threshold reso- 
nances {Su, -Pn, Pi 3} plus the sub-threshold Pi3(1720) 
resonance yields the best fit to all the available rj pro- 
duction data in the energy range considered in this work 
for photon- and hadron-induced reactions. The sub- 
threshold Pi3(1720) resonance, in particular, is required 
to help reproduce the observed shape of the rj angu- 
lar distribution at the excess energy of Q = 46.6 MeV 
measured by the COSY-11 Collaboration [Io|. We re- 
call here that in the present analysis, the mass of the 
sub-threshold P i3 as well as its width were fixed at the 
outset at the centroid values of the Pi 3 (1720)**** quoted 
in PDG [3]. As has been mentioned in Sec. [Til we a ls° 
note that the NN — > NNrj' reaction rules out the set of 
above-threshold resonances {5n, Pn, -D13} which fits the 
photoproduction data as well as the set {Su, Pn, P13}. 
As discussed later in connection with the invariant pp 
mass distribution in Fig. the set {Su, Pn, P13} is 
unable to reproduce the measured invariant pp mass dis- 
tribution. In the following, for simplicity, we restrict 
ourselves to the results obtained in conjunction with the 
CLAS photoproduction data (cf. Table |H column labeled 
"free p" CLAS) because the results corresponding to the 
other parameter sets shown in Table Q] can fit the NNrj' 
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TABLE III. The (hadronic) parameter values as determined from the combined fit to the 'yN — > n'N (parameter values given 
in Table HI free p CLAS), tvN -> rj'N, and NN —> NNrj' reaction data. The values in boldface were kept fixed during the fit 
procedure. The values in square brackets for the branching ratios are the PDG quotes, (fljvjvrj' , A) = (1.00,0.53). 



parameters 


Si i (1925) 


Pn(2130) 


Pi 3 (1720) 


Pis (2050) 


M R ( MeV) 


1924 


2129 


1720 


2050 


T R (MeV) 


112 


205 


200 


140 


/W ( % ) 


6 


3 


0.09 


2 


Pn, (%) 


22 


1 


[11 ± 3] 16 


4 


Pn v (%) 


25 


[61 ± 60] 1 


[4.0 ± 1.0] 


8 


Pn p (%) 


25 


77 


[70-85] 76 


59 


Pn„ (%) 


22 


18 


7 


26 







(0.68, 1.00) 


(1.77, 1.00) 


(1.20, -) 


(1.38, -) 


(SflJViri A) 




(-0.35, 0.97) 


(0.23, 0.34) 


(-0.17, -) 


(-0.05, -) 




(3) \ 

9rn p ) 


(-0.74, 0.90) 


(0.63, 0.52) 


(0.27, — ) 


(0.61, -) 


fa (1) o (2) 


(-2.69, 0.06, — ) 


(2.90, -0.10, — ) 


(-14.97, 43.19, 8.86) 


(2.52, 8.09, 28.39) 


fo (1) o (2) 


> 9rNuj) 


(1.93, -1.46, — ) 


(2.38, 1.25, — ) 


(-4.70, 8.23, -41.61) 


(-1.76, -16.29, -23.43) 



data equally well. 

The set of the parameter values of the present model 
that directly affect the hadronic processes is displayed in 
Table IIII1 The values result from the combined fit of the 
rjf production in photon- and hadron-induced reactions. 
The resonance partial decay widths were calculated by 
folding the partial decay widths for a given decaying reso- 
nance mass and a given emitted meson mass with the cor- 
responding mass distributions. For the latter, we assume 
Gaussian distributions with widths given by the corre- 
sponding total widths of the resonance and of the meson. 
The branching ratio /3nt]> for the sub-threshold Pi3(1720) 
resonance is extremely small because this arises only from 
the far (upper) tail of its mass distribution. As such, it 
is subject to considerable uncertainties. In the present 
work, the resonance coupling constants cannot be deter- 
mined uniquely because the available data for rj' pro- 
duction are not sufficient to impose more stringent con- 
straints. In particular, the lack of pn — > pnn' data — 
there exist only three upper-limit total cross-section data 
points (see Fig.|9|) — makes it difficult to constrain quan- 
titatively the relative contributions of the isoscalar (rj, ui) 
and isovector (ir, p) meson exchanges. It is also clear that 
one needs to consider meson-production reactions other 
than n' production to better constrain the relevant cou- 
pling constants. 

Our result for the pp —> ppn' total cross section is 
shown in Fig. [9] We see that the data are nicely re- 
produced over a wide range of excess energy. The dy- 
namical content of the present model is also displayed. 
One sees that the spin- 1/2 and -3/2 resonance contri- 
butions [in Fig. 1 N u = 5 n (1925) + Pn(2130) and 
Ni3 = Pi 3 (1720) + Pi 3 (2050), respectively] have differ- 
ent energy dependencies. In the lower excess-energy re- 
gion, the spin- 1/2 resonance contribution is only slightly 



smaller than the spin-3/2 resonance contribution, but as 
the energy increases, the spin-3/2 resonance contribu- 
tion starts to dominate. Here, the dominant spin-3/2 
resonance contribution is due to the Pi3(2050), while 
the dominant spin- 1/2 resonance contribution is from the 
Sii(1925) resonance. Although the overall relative Nu 
and N±3 resonance-set contributions are well determined, 
the individual spin- 1/2 and -3/2 resonance contributions 
within the set are not well constrained by the existing 
data, which is the reason why we only show the contribu- 
tions of the sums of the spin-1/2 (Nu) and the spin-3/2 
(Ni 3 ) resonances. The nucleonic and mesonic currents 
yield contributions that are much smaller than those of 
the resonances in the entire excess-energy range shown in 
Fig. HI The result for the pn — > pnn' total cross section 
is also shown in Fig. |H1 The pnrj 1 together with the ppn' 
reactions helps constraining the isoscalar- and isovector- 
meson couplings (here, M = tt, n, p, to) to the resonances. 
Unfortunately, only the upper limit of the cross section 
in a limited energy range is currently available for the 
pnr]' process. Therefore, the isoscalar-isovector meson 
exchange content of the present model is subject to this 
limitation in the existing data. 

The results for the n' angular distribution in pp — > ppn' 
at the excess energies of Q = 46.6 and 143.8 MeV are 
shown in Fig. 1101 The data are reproduced very well. 
We recall that the sub-threshold Pi3(1720) resonance is 
needed to reproduce the experimentally observed rf an- 
gular distribution, especially, at Q — 46.6 MeV. Here, 
one might regard this resonance as simulating some miss- 
ing background in the present model. In this connection, 
however, we mention that this resonance is the closest 
known resonance to threshold that helps to reproduce 
the measured angular distribution. In the left-hand-side 
column, at Q = 46.6 MeV, we see that the angular dis- 
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FIG. 9. (Color online) Total cross section for NN -> NNr)' 
as a function of the excess energy Q = y/s — ^/so, where 
y/so = 2mjv + m v i. The results (blue solid curves) corre- 
spond to the parameter set determined in conjunction with 
the fit to the 7iV — > r/'N and nN — >• rj 'N data. For ppr)', the 
individual current contributions are also show: nucleonic cur- 
rent (red dashed curve), mesonic current (green dash-dotted), 
TVii = 5*11(1925) + Pn(2130) resonance current (magenta 
dash-double-dotted), N 13 = Pi 3 (1720) + Pi 3 (2050) (maroon 
dotted). The ppr/' data are from Refs. [q-[l01|; the pnrj' data, 
which are upper limits, are from Ref. [llj. 



tribution is dominated by the 77' in the s-wave due to 
the completely dominant transition 3 Pq — > 1 Sqs. At 
the higher energy of Q = 143.8 MeV, the s-wave con- 
tribution still dominates to a large extent, leading to a 
nearly flat angular distribution as exhibited by the data, 
even though the 3 Pq — > 1 Sos partial- wave contribution is 
somewhat smaller at this high energy than at the lower 
energy of Q — 46.6 MeV. As can be seen in the panels 
on the right-hand-side column in Fig. [TU1 the flat angu- 
lar distribution, especially at higher energy, is achieved 
in the present model by an interference among different 
currents, in particular, between the spin-1/2 and -3/2 
resonance currents. 

In Fig. 111! the result for the prj' invariant mass distri- 
butions at an excess energy of Q = 16.4 MeV is shown. 
The data are reproduced well. 

The pp invariant mass distribution together with the rj 
angular distribution poses a relatively strict constraint on 
the set of resonances, provided they are above-threshold 
resonances. In particular, the set of above-threshold res- 
onances {Su, Pu, -D13}, which describes the photopro- 
duction data as well as the set {Sn,Pn, Pis}, is unable 
to reproduce the higher energy region of the measured 



pp invariant mass distribution by the COSY-11 collabo- 
ration [7J. 

As we shall discuss below, the present result for the 
pp invariant mass distribution in pp — > ppr]' reaction has 
an interesting implication on the issue of the reaction 
mechanisms in the pp — !• ppr/ reaction. In the latter reac- 
tion, there has been observed a significant enhancement 
of the cross section for larger pp invariant mass values 
compared to that given by the phase-space plus the pp 
FSI |20l . l2~j| . One of the possible explanations for this en- 
hancement is the relatively strong Nr) FSI. However, as 
pointed out in Ref. . ?| , the shape of the pp invariant mass 
distribution data in ppr\ at Q = 15.5 MeV is practically 
the same as that of pprj' shown in Fig. [T2] Since the Nr]' 
FSI is much smaller than the Nr) FSI, the explanation of 
the enhancement based on the Nr) FSI was ruled out in 
pp — > ppr) in Ref. Q- In Ref. [Hj], an alternative expla- 
nation, based on the higher partial-wave (final state P- 
wave) contribution, was proposed, together with a way to 
verify this proposed mechanism in a model-independent 
manner. Yet another explanation is the energy depen- 
dence of the basic production amplitude J [cf. Eq. (|A9I) ] 
as proposed in Ref. 23]. 

Figure [12] shows our results for the pp invariant mass 
distribution in pp — > ppr/ revealing a good agreement 
with the COSY-11 data 0. In the left panel of Fig. QJ 
one sees that the pp invariant mass distribution is prac- 
tically exhausted by the 3 Pq —> 1 S'os partial wave. This 
is quite surprising in view of the findings of Ref. [22| 
mentioned above, where a significant final-state P-wave 
contribution was found in the higher pp invariant mass 
region in the pp — > ppr) reaction. The present finding im- 
plies that the S'-wave basic production amplitude (J) in 
the present model should have an energy dependence as 
proposed in Ref. [23| , since the pp invariant-mass depen- 
dence introduced by the pp FSI is not enough to account 
for the enhancement of the measured pp invariant-mass 
distribution at larger invariant masses. This finding tells 
us that the conclusion reached in Ref. [7| ruling out the 
Nr) FSI as a possible source of the enhancement in the pp 
invariant mass distribution at larger invariant mass val- 
ues based on the comparison of the corresponding shapes 
in ppr) and ppr]' has to be taken with caution since there 
might be different mechanisms operating in these reac- 
tions as shown explicitly here in Fig. [T2] At this stage, 
it is natural to ask what the underlying dynamics is in 
the S'-wave contribution that accounts for the enhance- 
ment of the pp invariant mass distribution at larger in- 
variant mass values in pp — > ppr)' as compared to that 
in the pp — > ppr) reaction, where the enhancement arises 
from the 1 Sq — > 3 Pqs partial wave. In the right panel 
of Fig. 1121 we show the individual current contribution 
to the pp invariant mass distribution. We see that the 
enhancement at higher values of invariant mass is largely 
due to the constructive interference between the spin- 
1/2 (green dash-dotted curve) and the spin-3/2 (maroon 
dash-double-dotted curve) resonance contributions. The 
present model prediction for pp invariant-mass distribu- 
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FIG. 10. (Color online) rj angular distribution in pp — > pprj' in the center-of-momentum frame of the system for two excess 
energies of Q — 46.6 and 143.8 MeV. The results (blue solid curves) correspond to the parameter set determined in conjunction 
with the fit to the rj photoproduction. The panels in left hand-side column show the contributions from the partial waves 
3 Po — > 1 Sos (red dashed curve), sum of the s-wave (Z = 0) rj (green dash-dotted), s +p-waves (Z = 0,1) (maroon dash- 
double-dotted), and s +p + d -waves (Z — 0,1,2) (magenta dotted). The panels in the right-hand-side column show the 
nucleonic+mesonic currents (red dashed), the spin- 1/2 resonances (green dash-dotted), the spin-3/2 resonances (maroon dash- 
double-dotted), and the sum of the spin- 1/2 and -3/2 resonances (magenta dotted) contributions. The data are from the 
COSY-11 collaboration (46.6 MeV) [5j and from DISTO (143.8 MeV) @|. 



tion at a higher excess energy of Q = 46.6 MeV is shown 
in the left panel of Fig. [T3] Here one sees an onset of the 
1 Z?2 — > 3 P2S contribution. This result together with the 
result in Fig. reveal that the P-wave contribution in 
pp — > pprj' is much smaller than that in pp — > ppr\- For 
a close comparison with the pp invariant mass distribu- 
tion in the pp — » pprj reaction at higher excess energies, 
we show in the right panel of Fig. |13 | t he corresponding 
results at Q = 41 MeV from Ref. [22|. Here we see a 
striking difference between the model results for the two 
reactions. In pp —> pprj, there is a very large enhancement 
in the pp invariant mass distribution at higher invariant 
mass values due to the 1 Sq — > 3 Pqs partial wave con- 
tribution. By contrast, this partial wave contribution is 
minimal in the pp — > pprj reaction. It would be very in- 
teresting to have measurements of the pp invariant mass 
distribution at this excess energy in pp — > pprj to ver- 
ify the present model prediction. Obviously, these are 
model-dependent results. As such, it will be very inter- 



esting to verify them in a mo del- independent manner as 
pointed out in Ref. p3 [. 



VI. ttN r)'N 

Experimental data for the irN — > -q'N reaction are 
scarce. The only available data are the total cross sec- 
tions for ir~p — > rjn and n + n — > rjp [24l - l26l | which 
are subject to large uncertainties, as can be seen in 
Fig. Q31 Notwithstanding the fact that these data of- 
fer little constraints for the model parameters, they were 
included in the global fit and the corresponding fit re- 
sults are also shown in Fig. [TH We note here that 
within the present model the results for 7r~p — > rjn and 
7r + n — > rjp are identical. An interesting feature of the 
present model result is the double-bump structure caused 
by the #11(1925) and P u (2130) resonances (cf. Table 
IIIip . The former is just about 20 MeV above thresh- 
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FIG. 11. (Color online) The pr\ invariant mass distributions 
in pp — > pprj' at an excess energy of Q — 16.4 MeV. The 
blue solid curve corresponds to the full result determined in 
conjunction with the fit to the rj photoproduction data. The 
data are from Ref. 0j- Here, the theoretical results have been 
multiplied by an arbitrary normalization factor of 1.478 in 
order to facilitate the comparison of the shape with the data. 



old. In view of their large uncertainties, the currently 
existing data shown in Fig. [14] can indeed accommodate 
such a structure, however, clearly more accurate data are 
needed for a definitive answer. If experimentally corrob- 
orated, such a bump structure would rule out the sub- 
threshold resonance-dominance assumption of Ref. [27| . 
where Sii(1535) resonance dominance is assumed to de- 
scribe both the ttN — > rj'N and NN -» NNij' cross- 
sections data since it is not possible to generate any bump 
structure from sub-threshold resonances alone. 



VII. SUMMARY 

In the present work, we have revisited the theoreti- 
cal description of rf production in photon- and nucleon- 
induced reactions to take into account the recent addi- 
tions of accurate data to the corresponding world data 
base 0-0, EH- All of the currently available data in 
the resonance-energy region considered in this work are 
nicely reproduced within the present model in a combined 
analysis of the reactions 7 A — > Nrf, NN — > NNrf, and 
ttN —> rj'N. Considering only spin-1/2 and -3/2 reso- 
nances, we have found that the data are reproduced with 
a minimum of four resonances, i.e. Pi 3 (1720), S , n(1925), 
Pn(2130), and Pi 3 (2050). The Pi 3 (1720) is a four-star 
resonance listed in PDG [l4j], and the later three res- 
onances can be tentatively identified with the two-star 
5u(1895), one-star P u (2100), and one-star Pi 3 (2040) as 
listed in PDG [l4| ■ All three above-threshold resonances 
quoted above are essential for achieving the fit quality 
obtained in this work. Leaving out any one of them de- 
teriorates the fit quality considerably. The high-precision 



CLAS photoproduction data [H constrain the masses of 
the above-threshold resonances very well. 

In the free ^p — > rj'p reaction, there is a signifi- 
cant discrepancy between the most recent CLAS [J and 
CBELSA/TAPS @ data. Since currently there is no 
clear reason to discard one set in favor of the other, these 
two data sets lead to differences in the extracted res- 
onance parameters depending on which set is used for 
the analysis. The major difference is in the extracted 
coupling strength of the Pn(2130) resonance, where the 
CBELSA/TAPS data |] yield a much larger coupling 
than the CLAS data Q due to the larger cross sections 
exhibited by the CBELSA/TAPS data at higher energies. 

The quasi-free 7p — > rj'p and rj'n reactions 

have been also considered as a part of the combined 
analysis of the photon- and nucleon-induced reactions. 
The latter (quasi-free) reaction helps to constrain the set 
of above-threshold resonances. Overall, these reactions 
were reasonably well described by folding the cross sec- 
tions of the corresponding free processes with the Fermi 
distribution of the nucleon inside the deuteron. The ra- 
tio of the neutron to proton electromagnetic couplings 
for the considered resonances were extracted. Overall, 
the CBELSA/TAPS proton quasi- free data [1] coincide 
with the proton free data within their uncertainties for 
most of the angles and energies, although the quasi-free 
data are more in line with the free CLAS data [|| than 
with the free CBELSA/TAPS data [|| at higher energies, 
where the effect of the Fermi motion of the nucleon in- 
side the deuteron is expected to have faded out. This 
might be an indication of a possible problem with the 
CBELSA/TAPS proton free data @. At lower energies, 
down to W — 1.935 GeV, the lowest energy for which 
the proton free data exist, the data of Ref. Q show no 
sign of the Fermi-motion effect, while the present model 
calculation exhibits this effect up to W ~ 2 GeV. Further 
studies are necessary to identify the cause of this seeming 
discrepancy between the data and the model result. 

We found that the existing photon-induced reac- 
tions themselves can be described by a set of three 
above-threshold resonances. However, the photoproduc- 
tion data alone cannot constrain the set of resonances 
uniquely. In the present analysis, we also incorporate the 
pp invariant mass distribution data in the NN — > NNrj' 
reaction to obtain a restriction to a single set of above- 
threshold resonances. The data clearly require the above- 
threshold Pi3(2050) resonance. Furthermore, the rf an- 
gular distribution data in this reaction cannot be de- 
scribed adequately without the sub-threshold Pi 3 (1720) 
resonance. Clearly, these findings illustrates that meson 
productions in NN collisions can help impose constraints 
on the resonances. As pointed out in Ref. [Hj], spin ob- 
servables in rj photoproduction, in particular the beam 
asymmetry, are much more sensitive to the details of the 
model than are the cross sections. We also expect that 
the analyzing power in NN — > NNrj' is sensitive to the 
excitation mechanism of a given resonance, as is the case 
for the NN — > NNrj reaction (28|. Also, in view of the 
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FIG. 12. (Color online) The pp invariant mass distributions in pp — > ppr/' at an excess energy of Q — 16.4 MeV. The blue 
solid curves correspond to the full result determined in conjunction with the fit to the rj photoproduction data. The left panel 
shows the 3 Po — SqS (red dashed curve) partial wave contributions. The right panel shows the nucleonic+mesonic currents 
(red dashed), the spin- 1/2 resonances (green dash-dotted), the spin-3/2 resonances (maroon dash-double-dotted), and the sum 
of the spin-1/2 and -3/2 resonances (magenta dotted) contributions. The data are from Ref. 0]. Here, the theoretical results 
have been multiplied by an arbitrary normalization factor of 1.478 in order to facilitate the comparison of the shape with the 
data. 




FIG. 13. (Color online) Left panel: prediction for the pp invariant mass distribution in pp — > pprj' at an excess energy of 
Q = 46.6 MeV. The blue solid curve corresponds to the result determined in conjunction with the fit to the rj' photoproduction 
data. The 3 Po — > 1 5'os (red dashed curve) and the 1 !^ — > 3 P2S (green dash-dotted curve) partial wave contributions are also 
shown. Right panel: the model result (blue solid curve) of Ref. [22] for the pp invariant mass distribution in pp — > pprj at an 
excess energy of Q = 41 MeV. The two dominant partial-wave contributions, 3 Po — > 1 Sos (red dashed curve) and 1 So — > 3 Pos 
(green dash-dotted curve), are also shown. The data are from Ref. [2(J. 



contrasting results in pp — > pprj and pp — > pprj' reactions 
shown in Fig. II 31 it will be very interesting to measure 
the pp invariant mass distribution at an excess energy of 
Q ~ 45 MeV in the latter reaction to learn more about 
the possible production mechanism(s) in this reaction. 
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Appendix A: Formalism 

The formalism used in the present work is the same 
as that of Refs. [H, [l2|. For completeness, we provide 
here a brief description of this approach whose dynamical 
content is summarized by the graphs of Figs. IT51 and ITT1 



1. Free photoproduction 

For the rf photoproduction, we employ the tree graphs 
of Fig. [15] with form factors at the vertices to account 
for their hadronic structure. The gauge invariance of 
this production current is ensured by a phenomenological 
contact current that accounts for the effects of the final- 
state interaction current not taken into account explicitly 
[Hill- Following Refs. [H,|30j], the generalized contact 
current is chosen to be 




FIG. 15. Diagrams contributing to 7p — > rfp. Time proceeds 
from right to left. The intermediate baryon states are denoted 
N for the nucleon, and R for the Sn and Pn resonances. The 
intermediate mesons in the i-channel are p and ui. The ex- 
ternal legs are labeled by the four-momenta of the respective 
particles and the labels s, u, and t of the hadronic vertices cor- 
respond to the off-shell Mandelstam variables of the respective 
intermediate particles. The three diagrams in the lower part 
of the diagram are transverse individually; the three diagrams 
in the upper part are made gauge-invariant by an appropri- 
ate choice (see text) of the contact current depicted in the 
top-right diagram. The nucleonic current (nuc) referred to in 
the text corresponds to the top line of diagrams; the meson- 
exchange current (mec) and resonance current contributions 
correspond, respectively, to the leftmost diagram and the two 
diagrams on the right of the bottom line of diagrams. 



given by 



u — p 



p 



(A2) 

with k, p, and p' denoting the four-momenta of the in- 
coming photon, initial nucleon and final nucleon, respec- 
tively, and 



F = l-h(l-5J s )(l-6 u f u ), 



(A3) 



with indices x — s,u corresponding to the Mandelstam 
variables appropriate for the respective kinematic situ- 
ations depicted in Fig. [15] The factors 5 X are unity if 
the corresponding channel contributes to the reaction in 
question, and they are zero otherwise; f x denotes the ap- 
propriate form factor. The parameter h may be an arbi- 
trary (complex) function, h — h(s, u, t) which, in general, 
is subject to crossing-symmetry constraints. (However, 
in the application discussed in this work, we simply take 
h as a fit constant.) 

The interaction Lagrangians as well as the form fac- 
tors that provide the meson-nucleon-baryon and photon- 
nucleon-baryon vertices involved in the amplitudes M^, 
M£, Mf, and M£ are given in Appendix [B] 



M? = r NNr ,,(q)C», (Al) 

where r^jvy stands for the NNr]' vertex operator, with- 
out the corresponding isospin operator and the form fac- 
tor. The latter enter the auxiliary current C M that is 



2. Quasi-free photoproduction 

Following Refs. [1, EM IMj] , the quasi- free rj' photopro- 
duction processes are described within a spectator model 
by folding the cross sections for the corresponding free 
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processes with the momentum distribution of the nucleon 
inside the deuteron. 

In the laboratory frame, the deuteron is at rest and the 
spectator nucleon inside the deuteron is on its mass-shell 
and has three- momentum p s = — p N . Therefore, the 
energy En of the participant nucleon inside the deuteron 
is given by 

En = M d - ^m 2 s +p 2 N , (A4) 

where m s denotes the mass of the spectator nucleon. 

The invariant mass square of the j(k) + N(pn) — > 
r]'(q) + N(p' N ) subsystem, Q 2 = (k + pn) 2 , can be ex- 
pressed in the laboratory frame, where the participant 
nucleon inside the deuteron has three- momentum p N , as 

Q 2 (p N ) = E%-p% + 2£ 7 (e n - p N ■ k) , (A5) 
I 



where E~ = (s — AfJ)/(2y / s) denotes the incident photon 
energy with s denoting the invariant mass squared of the 
7c? system and, Md, the deuteron mass, k = k/|k| with 
k denoting the three-momentum of the incident photon. 

Of course, one must have Q{p N ) > m N + m v > for 
the quasi-free ^N — > r/'N process to take place. To- 
gether with this condition, the invariant mass squared, 
Q 2 , in Eq. (|A5|) can be expressed in terms of the four- 
momentum transfer square t = {pd — p s ) 2 — p% = 
E N — p N , where pd and p s stand for the four-momenta 
of the deuteron and the spectator nucleon, respectively, 
as 

Q 2 {p N ) =t + 2E y ^p 2 N +t - p N ■ kj 

= Q 2 (t,p N ). (A6) 

The differential cross section of the quasi-free photo- 
production process is, then, approximated as 



da 
dQ 



(W,0) 



quasi 



d 3 p N \^(p N )\ 2 Q(Q{t,p N ) -m N -77ty) ^ 



(W' = Q(t,p 



N Ji 



(A7) 



free 



where W = Q(t = m 2 N , p N = 0). Q(x) is the usual step function that is equal to 1 for x > and otherwise. It has 
been introduce for the sole purpose of making explicit that Q(t, p N ) > rn n + m n for the free 7 + N — » r/'N process 
to take place. 6 is the scattering angle between outgoing meson and incoming photon, ^(pjy) is the deuteron wave 
function in momentum space. d<T/df2|f roc is the differential cross section for a free photoproduction process. 

A variant of Eq. (|A8[) is to restrict the participant nucleon to be on its mass-shell in the argument of dcr/dO|f ree in 
Eq. (JA7J) HEEa, i.e., 



da 

dJi 



(W,9) 



quasi 



d 3 PJ v|*(Piv)l 6(Q(i,p 



N I 



m N 



v) 



da 

dn 



(W' = Q(m 2 N , p N ), 



(A8) 



free 



There are at least two reasons for this on-shell restriction. 
One is the fact that the free cross sections that enter 
the equations above are "on shell". The other is that 
by restricting t — m^, the effect of the Fermi folding is 
to smear out the free cross section keeping the centroid 
position, corresponding to da /dQ\ hcc (W, 0), not to be 
shifted. The latter feature seems to better reproduce the 
quasi- free i] photoproduction data [l~8|, [l9[ . 

In the present work we employ Eq. (|A8|) . but also ob- 
tained results using Eq. (|A7[) and we find little difference 
between the two prescriptions. 



3. Hadronic production 

The ttN — > r/'N reaction is described here within a 
tree-level approximation, analogous to the description of 
ttN — > rjN in Ref. [l6j . We take into account the nucle- 
onic and resonance contributions as depicted in Fig. 1161 
In principle, one could also include the ^-channel dia- 
grams such as the (rank-two) tensor meson ai (1320) ex- 



change, whose decay branching ratio to ri'n is quoted to 
be BR{a 2 -> tj'tt) = 5.3 ± 0.9 x 10~ 3 >4]. We have 
not considered such contributions in the present work 
since we do not expect that including them would alter 
our results in any significant manner. The propagators, 
vertices, and form factors necessary for calculating the 
Feynman diagrams in Fig. \W\ are given in Appendix [B] 

The hadronic reaction NN — > NNrj' is described ac- 
cording to the model put forward in Refs. [16|, l28|, |32j. 
The DWBA amplitude M for this process is given by 

M = (1 + T f G f )J(l + GiTi), (A9) 

where T„, with n = i, f, denotes the NN T-matrix in- 
teraction in the initial (i) or final (/) state, and G n 
is the corresponding two-nucleon propagator (which ab- 
sorbs the factor i found in the DWBA formula given in 
Ref. [28j). J sums up the basic rj production mechanisms 
depicted in Fig. [TTJ The interaction Lagrangians as well 
as the form factors necessary for constructing the basic 
production amplitude J are given in Appendix B. In the 
absence of models capable of providing a reliable off-shell 
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FIG. 16. Feynman diagrams contributing to ttN — > ry'TV. The 
notation is the same as in Fig. 1151 
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NN initial state interaction (ISI), we consider it only 
in the on-shell approximation following Refs. [HI l33j . 
This was shown to be a reasonable approximation for 
calculating cross sections [Hj]. For the on-shell NN in- 
teraction, we consider the phase-shifts and inelasticities 
from the SAID partial- wave analysis [34{ . All the partial 
waves up to total angular momentum J — 7 are consid- 
ered. The NN FSI is generated using the Paris potential 
1351 where the Coulomb interaction is taken into account 



fully as described in Ref. [32J. We also follow Ref. [36| 
to convert the Paris NN interaction, which obeys the 
non-relativistic Lippman-Schwinger equation, to the one 
obeying the three dimensionally reduced version (a la 
Blankenbecler-Sugar) of the relativistic Bethe-Salpeter 
equation in order to be consistent with the relativistic 
covariant approach used in the present work. We also use 
the Blankenbecler-Sugar propagator for the two-nucleon 
propagator Gf in Eq. (|A9j) for consistency. 



Appendix B: Lagrangians, Form Factors, and 
Propagators 

The interaction Lagrangian used to construct our 
model for the basic production amplitudes is given be- 
low. For further convenience, we define the operators 



■■(+) 



75 



and 



r(-) = i. 



(Bl) 



1. Electromagnetic interaction Lagrangians 



c 



r , 



eN 



K 



-o» v d v A 



2M N 

e^e a ^ (d a A») (0Y) pi, 



9n'< 



N, (B2) 
(B3) 
(B4) 



where e stands for the elementary charge unit, and 
e = (l + r 3 )/2 and k = K p (l + r 3 )/2 + K n (l- r 3 )/2, with 



FIG. 17. Basic production mechanisms for AW — > NNrj' . 
Time proceeds from right to left. The full amplitude, with 
additional initial- and final-state contributions, is given by 
Eq. (|A9fl . As in Fig. 1151 N and R denote the intermediate 
nucleon and resonances, respectively, and M incorporates all 
exchanges of mesons n, r\, p, uj, a, and ao (= 8) for the nucleon 
graphs and n, p, and ui for the resonance graphs. External 
legs are labeled by the four-momenta of the respective par- 
ticles; the hadronic vertices s, u, and t here correspond to 
the same kinematic situations, respectively, as those identi- 
fied similarly in Fig. 1151 The nucleonic (nuc), resonance, and 
meson- exchange (mec) contributions referred to in the text 
correspond, respectively, to the first, second, and third lines 
of the diagrams on the right-hand side. 



the anomalous magnetic moments k p — 1.793 of the pro- 
ton and K n = —1.913 of the neutron; Mm stands for the 
nucleon mass. F M „ = d^A v — d v A^ with A^ denoting the 
electromagnetic field and e all \ v is the totally antisym- 
metric Levi-Civita tensor with e 0123 = +1. The meson- 
meson electromagnetic transition coupling constants in 
the above Lagrangians, g^'fr/ — 1-25 and 3j/"7 — 

0.44, 

are extracted from a systematic analysis of the radiative 
decay of pseudoscalar and vector mesons based on fla- 
vor SU(3) symmetry considerations in conjunction with 
vector-meson dominance arguments [32^ : their signs are 
inferred in conjunction with the sign of the coupling con- 
stant <7 7r „ 7 (v — p, lj) determined from a study of pion 
photoproduction in the 1 GeV energy region [37]. The 
resulting rj'vj vertex is multiplied by the form factor 
f v (t) which describes the off-shell behavior of the inter- 
mediate vector meson with squared momentum transfer 
t = (p — p') 2 (cf. fourth diagram in Fig. [15]). We use the 
dipole form 



fv(t) 



A* 



A* 



(B5) 



The cutoff A*, taken to be identical for both p and lu, is 
a fit parameter. 

The resonance-nucleon-photon-transition Lagrangians 
are 



o (1) 
'~2Mn 



RT^(r^ v (d v A^)N + H.c, (B6a) 
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-RNj 



o (1) 



2M 

(2) 



N 



R^T {±) F^d v N + H.c. 



(B6b) 



where the superscript of £rn~/ denotes the spin and par- 
ity of the resonance R. The coupling constants g RNl 
(i = 1, 2) are fit parameters. 



C-nnv = —9nnvN 



2M 



N 



V» } N, (B7c) 



where the parameter A was introduced in Cnnp to in- 
terpolate between the pseudovector (A = 0) and the 
pseudoscalar (A = 1) couplings. The NNrj' coupling 
constant, gNNrj'i * s a ^ parameter. All the other cou- 
pling constants in the above Lagrangians are taken from 
Ref. [38f| . with the exception of gNNu — 10 as explained 
in RefiUl. 

The rj'vv (v = p, ui) Lagrangians are 



2. Hadronic interaction Lagrangians 



In this subsection, we use S (= c, <2o), P (= rj 1 if), and 
Vfi (= Wp,^) to denote the scalar, pseudoscalar, and 
vector meson fields, respectively. The vector notation 
refers to the isospin space. For isovector mesons, S = 



P ■ t, and Vu 



The Lagrangians for meson-nucleon interactions are 



£nns = 9nns NNS, 
Cnnp = -g NNP N \ 



iX + 



1 - A 

2M 



N 



(B7a) 
P}N, (B7b) 



9-n'pp 



r 



9v'< 



2m 



Tj' 



^ (d a pn (flY) ■ p\ (B8) 
z^xu (d a ^) (d x r)') u". (B9) 



The coupling constants, g^i pp = 5.51 and g v > uu — 5.42, 
are obtained from a systematic analysis of the radiative 
decay of pseudoscalar and vector mesons based on SU (3) 
symmetry considerations; their signs are inferred, in 
conjunction with vector-meson dominance assumptions, 
from the sign of the coupling constant g nvi (v = p, uj) 
determined from a study of pion photoproduction in the 
1 GeV energy region [37j . 



For nucleon resonances, 



C 



RNP 



C 



C 



ah 

RNV 

(§*) 
RNP 



C 



(§*) 
RNV 



-9rnp r \ r 



±iX 



1 - A 



M R ± M 



A' 



gpNv 



2M 



N 



RT^a^ {d v V»)N + H.c. 



^ R f^R it T^(d"P)N + B.c, 
Alp 



o (1) 

9 RN V 

2M. 



,(2) 



P^iV + H.c, 



(3) 



n 4:M N 4M N 



(BlOa) 
(BlOb) 
(BlOc) 

(BlOd) 



where the Lagrangian (|B10a[) contains the pseudoscalar-pseudovector mixing parameter A, similar to Eq. (IB7bl) . 

I 



Each hadronic vertex obtained from the interaction La- 
grangians given in this subsection is multiplied by a phe- 
nomenological cutoff function 



f(p' 2 ,p 2 ,q 2 ) = fB(p' 2 )fB(p 2 )fM(q 2 ), 



(Bll) 



with n = 1 for a scalar or a pseudoscalar meson and n = 2 
for a vector meson, tum denotes the mass of meson M . 
The values of Am are taken the same as those used in 
Ref. [H. 



where p' and p denote the four-momenta of the two 
baryons, and q is the four-momentum of the meson at 
the three-point vertex. Here, we use 



f B (x) 



A% + (x- Mir 



(B12) 



where the cutoff Kb — 1200 MeV is taken the same for 
all the baryons B, and /m(<Z 2 ) is given by 



fM(q 2 ) 



A 2 



'■M 



A 2 



(B13) 



3. Energy-dependent resonance widths 

Our formalism is adapted to accommodate energy- 
dependent resonance widths with the appropriate thresh- 
old behavior. 

For a spin- 1/2 resonance propagator, we use the ansatz 



1 



m R 



m R 



tL 



P 



(B14) 

where uir is the mass of the resonance with four- 
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momentum p. T is the width function whose functional 
behavior will be given below. 

For spin-3/2, the resonant propagator reads in a 
schematic matrix notation 



S 3 /2(P) 



{i>-m R )g-i—T 



A. 



(B15) 



All indices are suppressed here, i.e., g is the metric ten- 
sor and A is the Rarita-Schwinger tensor written in full 
detail as 



AF - 
— 



-g 



i 



Lt v 



2p^p v 



S/3a + 



P' 



3 ,p£ ,tu ■ 3m R ■ 3m R 

(B16) 

where /3, a, and e enumerate the four indices of the 7- 
matrix components (summation over e is implied). The 
inversion in (|B15I) is to be understood on the full 16- 
dimensional space of the four Lorentz indices and the 
four components of the gamma matrices. 



In both cases, we write the width T as a function of 
W = y/s according to 



r(w) = r R J2f3 l r l (w) 7 



(B17) 



where the sums over i accounts for decays of the reso- 
nance into two- or three-hadron channels and into radia- 
tive decay channels. The total static resonance width is 
denoted by T R and the numerical factors /3j describes the 
branching ratios into the various decay channels, i.e., 



JV 



i=l 



(B18) 



Similar to Refs. [39|-|42|, we parameterize the width func- 
tions Ti (which is normalized to unity at W — m R ) to 
provide the correct respective threshold behaviors. The 
details may be found in Ref. fl2l |. 
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